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Novel phase of beryllium fluoride at high pressure 


Maksim S. Rakitin/’[^ Artem R. Oganov,^’Haiyang M. Mahdi Davari 

Esfahani/ Xiang-Feng Zhou/ Guang-Rui Qian/ and Vladimir L. Solozhenko^ 

^Department of Geoseienees, State University of New York, Stony Brook, NY 11794, USA 
^Skolkovo Institute of Seienee and Teehnology, Skolkovo Innovation Center, Bldg. 3, Moseow 143026, Russia 
^Moseow Institution of Physies and Teehnology, 9 Institutskiy Lane, 

Dolgoprudny City, Moseow Region If 1700, Russia 
"^Sehool of Materials Seienee, Northwestern Polyteehnieal University, Xi’an 710072, China 
^Shenyang National Laboratory for Materials Seienee, Institute of Metal Researeh, 

Chinese Aeademy of Seienees, Shenyang 110016, China 
^LSPM-CNRS, Universite Paris Nord, 93430 Villetaneuse, Franee 
(Dated: June 15, 2015) 

A previously unknown thermodynamically stable high-pressure phase of BeF 2 has been predicted 
using the evolutionary algorithm USPEX. This phase occurs in the pressure range 18-27 GPa. Its 
structure has C2/c space group symmetry and contains 18 atoms in the primitive unit cell. Given 
the analogy between BeF 2 and Si 02 , silica phases have been investigated as well, but the new phase 
has not been observed to be thermodynamically stable for this system. However, it is found to be 
metastable and to have comparable energy to the known metastable phases of Si 02 , suggesting a 
possibility of its synthesis. 

PACS numbers: 


I. INTRODUCTION 

Beryllium fluoride has many applications, such as 
coolant component in molten salt nuclear reactor^^^, 
production of special glasse^^, manufacture of pure 
berylliunP, etc. Structurally, BeF2 phases are similar 
to the phases of Si02 (Fig. [^: a-quartz phase of BeF2 
and Si02 is stable from 0 to ^2 GPa, and then trans¬ 
forms to coesite phase which persists up to ^8 GPa, 
and then transforms to stishovite (rutile-type phase) in 
SiOg!^. However, the behavior of BeF2 experimentally is 
not known for pressures above 8 GPa (see Scheme 1 in 
RefPl). 
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FIG. 1. Phase diagrams of Si02P and BePj^at low (up to 
room) temperatures. 

One of our goals in present paper is to reveal which 
phase transitions can occur at higher pressures in BeF2. 
Beryllium compounds are extremely toxic for humans, 
and this limits experimentation. Gomputer simulation 
is a safe and cheap alternative to investigate such struc¬ 
tures. In recent ah initio stud}!^ authors explored 13 
well-known AB2 structure types for their possible sta¬ 
bility for BeF2: a-quartz-type (P3i21), /-quartz-type 
(P6222), Qf-cristobalite-type (P4i2i2), /-cristobalite- 
type (Pd-3m), cubic GaF2-type (Pm-3m), Qf-PbGl2-type 
{Pnma), Ni2ln-type {P6s/mmc), coesite-type {C^jc), 
rutile-type (P42/mnm), baddeleyite-type (P2i/c), a- 


Pb02-type {Pbcn), Qf-GaGl2-type {Pnnm) and pyrite- 
type {Pa-3) structures. They found that the sequence 
of pressure-induced phase transitions of BeF2 up to 
50 GPa is as follows: a-quartz-type GPa^ coesite- 

type-rutile-type-^ Qf-Pb02-type struc¬ 

tures. Although BeF2 under pressure has been theoret¬ 
ically investigated by Yu et alW, we revisit these results 
to check for previously unknown structure(s), and we ex¬ 
plore the relevance of these findings for Si02. 


II. COMPUTATIONAL DETAILS 

Gomputer simulations of BeF2 and Si02 has been per¬ 
formed in two steps: (1) prediction of a new structure 
of BeF2 using USPEX evolutionary algorithm; (2) calcu¬ 
lation of properties of BeF2 and Si02 in the wide range 
of pressures from 0 to 50 GPa with a 1 GPa step using 
DFT. 

To find stable lowest-energy crystals structures, we 
performed fixed-composition search of the BeF2 system 
at different pressures (15, 20 and 25 GPa) using the US¬ 
PEX code^Wm, in conjunction with first-principles struc¬ 
ture relaxations using density functional theory (DFT) 
within the Perdew-Burke-Ernzerhof (PBE) generalized 
gradient approximation (GGA)P^, as implemented in the 
VASP packagj^. We employed projector augmented 
wave (PAW)P^ potentials with 2 valence electrons for Be 
and 7 — for F. The wave functions were expanded in 
a plane-wave basis set with the kinetic energy cutoff of 
600 eV and F-centered meshes for Brillouin zone sam¬ 
pling with reciprocal space resolution of 27rx0. 10 A-k 

We used the VASP package to carefully reoptimize 
the obtained structures before calculating phonons, elas¬ 
ticity, electronic density of states (DOS), hardness of 
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BeF2 and Si02. For these relaxations, we also used the 
plane-wave cutoff of 600 eV and k-meshes with resolu¬ 
tion of 0.10 Phonons calculations have been per¬ 

formed using Phonop^nSl and Quantum EspressJ^ codes 
for the relaxed structures at pressures where these struc¬ 
tures are found to be thermodynamically stable. Hard¬ 
ness was calculated using 3 methods: Lyakhov-Oganov 
modeP^ based on the strength of bonds between atoms 
and bond network topology, Chen-Niu modeP^ which 
uses elastic constants obtained from DFT calculations 
and Mukhanov-Kurakevych-Solozhenko thermodynamic 
model of hardnes^^. 


III. RESULTS AND DISCUSSION 

USPEX allowed us to find a new structure of BeE2, 
stable at 18-27 GPa (Eig. |^. The structure has C2/c 
space group and contains 12 formula units in the Bravais 
cell (6 in the primitive cell) with a=8.742 A, 6=8.695 A, 
c=4.178 A and /3=66.1° (at 20 GPa). Galculated den¬ 
sity of this new C2/c phase is 4.2% higher than density 
of coesite phase, both at 20 GPa. Eor reference, here 
are lattice parameters for BeE2-stishovite at 30 GPa: 
a=6=3.986 A, c=2.501 A and a=/3=j=90°. The value 
of the bulk modulus Ho=22.4 GPa of the C2/c structure 
of BeE2 with its pressure derivative Hq=3.9 was obtained 
from a least-snares fit using the Murnaghan equation of 
stat^^ (Fig. 1^. The zero-pressure unit cell volume was 
taken as Vb=213.7 A^. 


FIG. 2. C2/c structure of BeF 2 , stable at 18-27 GPa. 


A. Thermodynamic stability 

We have calculated the enthalpies of a-quartz (P3221), 
coesite ((72/c), coesite-II ((72/c), stishovite (P42/mnm), 
a-Pb02-type {Pbcn) structure and our new structure 
((72/c) for both BeF2 and Si02 at different pressures 
from 0 to 50 GPa with a 1 GPa step. The results are 
presented in Fig. 



FIG. 3. Equation of state of BeF 2 C2/c structure. 


1. BeF 2 underpressure 

Eor the case of BeE2 a-quartz structure is stable from 
0 to 4 GPa, followed by coesite structure stable from 4 to 
18 GPa, and the (72/c structure is found to be stable be¬ 
tween 18 and 27 GPa, which then gi ves p lace to stishovite 
structure at higher pressures (Eig. |4(a)] ). We see transi¬ 
tion from coesite-type to (72/c, then to rutile-type, but at 
much higher pressure (27 GPa against 6.47 GPa in Ref.l^, 
where LDA was used). Accordin^o Demuth et the 
LDA approximation used in Ref.l^ underestimates phase 
transition pressures, whereas using the GGA yields more 
reliable results. The (a-Pb02-type structure is not sta¬ 
ble at any pressure (in the investigated interval from 0 
to 50 GPa) for BeE2 (though it is close to stability at 
^27 GPa), while for Si02 it is indeed stable at pressures 
above ~80-90 GPeP. 


2. Si 02 under pressure 

Erom Eig. |4(b)] it is clearly seen that in Si02 the tran¬ 
sition from a-quartz to coesite occurs at 5 GPa, followed 
by transformation to stishovite at ^7 GPa, which contin¬ 
ues to be stable up to 50 GPa. This phase transition se¬ 
quence is in good agreement with experiment^ and with 
the GGA results by Demuth et Oganov et al^ and 
LDA results of Tsuchiya et it is known thouglP^ 

that the GGA is more accurate than the LDA for phase 
transition pressures. The new structure is not stable at 
any pressure for Si02, but at 0 GPa is only 3.4 meV/atom 
higher in energy than a-quartz, and should be synthesiz- 
able as a metastable phase. Our results of coesite ^ 
coesite-II transition are in good agreement with recent 
study of Cernok et where they observe coesite at 

20.3 GPa, and after an abrupt change in the diffraction 
pattern between ^20 and ^28 GPa — coesite-II at 27.5 
and 30.9 GPa. 
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FIG. 4. Enthalpies (relative to a-quartz) Qf |(a)| BeF 2 and |(b)| 
Si 02 phases as a function of pressure. 


S. Metastahle structures of Si02 


It is well known that Si02 n-quartz is thermodynami¬ 
cally stable at ambient pressure. However, there are nu¬ 
merous known Si 02 polymorphs which are metastable, 
but exist in nature or can be synthesized. We examined 
Si 02 feldspar, baddeleyite, melanophlogite and mogan- 
ite at 0 GPa. El Goresy et al^ claimed a baddeleyite- 
like post-stishovite phase of silica in the Shergotty me¬ 
teorite, however later that controversial phase turned 
out to be (a-Pb 02 -like silicgP^. Our calculations confirm 
that the baddeleyite-like form of Si 02 is very unfavor¬ 
able at 0 GPa and spontaneously (barrierlessly) trans¬ 
forms into the (a-Pb02-like structure. We have found 
that Si02-feldspar, moganite and melanophlogite are en¬ 
ergetically very close to the stable phase (a-quartz) and 
to the new C2/c structure. Differences in enthalpy be¬ 
tween melanophlogite, new structure and a-quar tz ar e 
less than 0.01 eV/atom (see zoomed inset in Fig. |4(b)] ). 
The fact that complex open structure of melanophlog¬ 
ite (138 atoms/cell) has a slightly lower energy than a- 


quartz, can be explained by errors of the GGA, which 
were discussed in details b y Dem uth et alW^. They also 
found /3-cristobalite (Fig. |4(b)[ ) is lower in energy by 
about 0.03 eV/Si02 than a-quartz, confirmed by calcu¬ 
lations of Zhang et showing that the GGA slightly 
overstabilizes low-density structures. 


B. Lattice dynamics 

Since the new structure of BeF 2 appears to be ther¬ 
modynamically stable, analysis of dynamical stability 
(phonon dispersion) has been performed for this struc¬ 
ture as well as for all other structures at pressures where 
they were found to be thermodynamically stable. Our 
results show that BeF 2 a-quartz at 0 GPa, coesite at 
5 GPa, new structure at 25 GPa and stishovite at 30 GPa 
do not have imaginary frequencies. Similar results are ob¬ 
served for Si02 a-quartz at 0 GPa, coesite at 5 GPa and 
stishovite at 10 GPa. Fig. shows dynamical stability of 
the new structure of BeF 2 since no imaginary frequencies 
are observed in the phonon dispersion plot. 



FIG. 5. Phonons dispersion curves showing dynamical stabil¬ 
ity of the C2/c structure of BeF 2 at 25 GPa. 


C. Electronic properties 

According to Fig. all BeF 2 phases are insulators, 
the DFT band gap increases from ^7 to ^10 eV with 
increasing pressure from 0 to 30 GPa and the value of 
the gap is in good agreement with data of Yu et al^. 

For Si02 (Fig. we also observe insulating behavior, 
and the band gap is about 6 eV and remains almost un¬ 
changed with increasing pressure. 


D. Hardness 

Three models have been exploited to calculate hard¬ 
nesses — the Lyakhov- Ogano\ff3, Chen-NirP^ and 
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FIG. 6. Density of states of BeF 2 in the a-quartz (at 0 GPa), 
coesite (at 5 GPa), C2/c structure (at 25 GPa), and stishovite 
(at 30 GPa) phases. 


Mukhanov-Kurakevych-SolozhenkcP^ models. First ap¬ 
proach is based on concepts of bond strengths and bond 
topology to compute hardness. Detailed description of 
the methodology can be found in RefP. This model 
has been implemented in the USPEX code, and for 
greater convenience has also been implemented as an 
online utility available at http: //han. ess . sunysb. edu/ 
hardness/, The second method of hardness calculation 
is Chen-Niu model, which is based on elastic tensor com¬ 
ponents and also implemented in the USPEX code. The 
third one is a thermodynamic model of hardness. 

The results can be seen in Table [l| Experimental data 
are provided where available — Vickers hardness of Si02- 
quartJ^, Si02-coesit^I^and Si02-stishovit^^. Erom Ta¬ 
ble [I] it is clearly seen that the calculated hardness of 
Si 02 quartz and stishovite is much higher than one of 
BeE 2 analogs. The hardness of BeE 2 and Si02 in the 
new C2/c structure is comparable with the hardness of 
Qf-quartz and coesite. 


IV. CONCLUSIONS 




(b) 



We have examined thermodynamic, vibrational, elec¬ 
tronic and elastic properties of BeE 2 and Si 02 phases us¬ 
ing DET calculations. The sequence of pressure-induced 
phase transitions of BeE 2 up to 50 GPa is as follows: a- 

quartz-type -coesite-type -^ G2/c -^ 

stishovite (rutile-type) structures. We found a new phase 
of BeE 2 which is thermodynamically stable at pressures 
from 18 to 27 GPa. This phase is not observed in Si02, 
but could be synthesized in principle. Electronic proper¬ 
ties analysis has shown BeE 2 and Si 02 remain insulating 
in a wide range of pressures (from 0 to 50 GPa). Hard¬ 
ness of BeE 2 and Si 02 in the new structure is comparable 
with hardness of a-quartz and coesite at 0 GPa. Hard¬ 
nesses of metastable Si 02 structures have been examined 
as well. 
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FIG. 7. De nsity of states of Si0 2 in the (a) o-quartz (at 
0 GPa), [(i^ coesite (at 5 GPa) and |(c)| stishovite (at 10 GPa) 
phases. 
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TABLE L Hardness of BeF 2 and Si02 structures at 0 GPa in GPa. For the metastable Si02 structures we present enthalpies 
relative to a-quartz (in eV per formula unit). 



Lyakhov- Oganov 

BeF2 

Ghen-Niu 

Mukhanov et al^ 

Lyakhov- 0 ganov 

Si02 

Chen-Niu 

Experiment 

Quartz 

7.1 

7.5 

11.0 

20.0 

12.5 


12.(1: 


Goesite 

8.2 

8.3 

11.7 

22.3 

8.4 


20.f 


New structure 

7.3 

6.8 

13.5 

19.1 

6.7 


— 


Stishovite 

8.2 

12.7 

15.1 

29.0 

28.7 


33.(0 



Metastable structures (Si02 only): 






Relative enthalpy, eV/f.u. 


Hardness, GPa 







Lyakhov-Oganov model 

Ghen-Niu model 


Feldspar 

0.047 

6.7 

11.8 



Baddeleyite 

0.726 

29.6 

28.0 



Melanophlogite 

-0.013 

12.5 

3.3 



Moganite 

0.003 

19.5 

12.8 




^ Thermodynamic model of hardness (Ref.^J) 
^ Vickers hardness 


ACKNOWLEDGMENTS 

We thank the National Science Foundation 
(EAR-1114313, DMR-1231586), DARPA (Grant 

No. W31P4Q1210008), the Government of Russian 
Federation (grant No. 14.A12.31.0003), and Foreign 
Talents Introduction and Academic Exchange Program 
(No. B08040). Also, we thank Dr. V.A. Mukhanov for 
valuable comments. 


Appendix A: Densities of BeF 2 and Si02 structures 

Table [TT] shows densities of BeF 2 structures at 0 and 
20 GPa and Si02 structures at 0 GPa. 


TABLE 11. Densities of BeF 2 and Si02 structures. 


System 

Number of 

atoms 

Volume, 

A/cell 

Density, 

g/cm'^ 

BeF 2 at 0 GPa: 

a-quartz 

9 

105.167 

2.244 

coesite 

24 

254.636 

2.472 

coesite-H 

96 

1021.960 

2.464 

C2/c 

18 

213.696 

2.209 

stishovite 

6 

47.771 

3.294 

BeF 2 at 20 GPa: 

a-quartz 

9 

73.078 

3.230 

coesite 

24 

202.001 

3.116 

C2lc 

18 

145.159 

3.252 

stishovite 

6 

41.492 

3.793 

Si02 at 0 GPa: 

a-quartz 

9 

116.934 

2.580 

coesite 

24 

283.341 

2.839 

coesite-H 

96 

1137.296 

2.830 

C2/c 

18 

243.569 

2.477 

stishovite 

6 

48.185 

4.174 

Q:-Pb02-type 

12 

94.623 

4.251 


Appendix B: GIF file of BeF 2 C2/c structure at 
20 GPa 


# GIF file 

# This file was generated by FINDSYM (H.T. Stokes) 
data_findsym-output 

_syininetry_space_group_nanie_H-M 1 2/c 
_syinmetry_Int_Tables_number 15 

_cell_length_a 8.74241 

_cell_length_b 8.69478 

_cell_length_c 4.17800 

_cell_angle_alpha 90.00000 

_cell_angle_beta 66.07927 

_cell_angle_gainma 90.00000 

loop_ 

_space_group_symop_operation_xyz 

x,y,z 

-x,y,-z+l/2 

-x,-y,-z 

x,-y,z+l/2 

x+l/2,y+l/2,z 

-x+l/2,y+l/2,-z+l/2 

-x+1/2,-y+1/2,-z 

x+l/2,-y+l/2,z+l/2 


loop_ 

_atom_site_label 

_atom_site_type_symbol 

_atom_site_fract_x 

_atom_site_fract_y 

_atom_site_fract_z 

_atom_site_occupancy 


Bel 

Be 

0.30175 

0.08755 

0.27471 

1.00000 

CN 

CD 

PQ 

Be 

0.00000 

0.18404 

0.25000 

1.00000 

FI 

F 

-0.11350 

0.09592 

0.11433 

1.00000 

F2 

F 

0.14707 

0.43232 

0.42082 

1.00000 

F3 

F 

-0.11681 

0.27080 

-0.42582 

1.00000 
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